Bose-Fermi mixtures of self-assembled filaments of fermionic polar molecules 
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Fermionic polar molecules in deep ID optical lattices may form self-assembled filaments when 
the electric dipoles are oriented along the lattice axis. These composites are bosons or fermions 
depending on the number of molecules per chain, leading to a peculiar and complex Bose-Fermi 
mixture, which we discuss in detail for the simplest case of a three-well potential. We show that 
the interplay between filament binding energy, transverse filament modes, and trimer Fermi energy 
leads to a rich variety of possible scenarios ranging from a degenerate Fermi gas of trimers to a 
binary mixture of two different types of bosonic dimers. We study the intriguing zero temperature 
and finite temperature physics of these composites for the particular case of an ideal filament gas 
loaded in fD sites, and discuss possible methods to probe these chain mixtures. 
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PACS numbers: 03.75.Lm,05.30. Jp 

A new generation of experiments is starting to explore 
systems where the dipole-dipole interaction (DDI) plays a 
significant and possibly dominant role. Due to the long- 
range anisotropic character of the DDI, dipolar quan- 
tum gases offer a fascinating novel physics [l], 0). Ex- 
citing phenomena have been recently reported for exper- 
iments on Bose-Einstein condensates (BECs) of atomic 
magnetic dipoles, especially on Chromium BEC I3|, but 
also on spinor Rubidium BECs 0], Potassium [5| and 
Lithium [6| . Magnetic atomic dipoles are however rather 
weak. On the contrary heteronuclear molecules, espe- 
cially at their lowest rovibrational state, may present a 
very large electric dipolc moment (> I Debye) @, d, HJ. 
Although quantum degeneracy has not been yet achieved, 
the rapid pace of development allows to expect degener- 
ate gases of polar molecules in the next future. These 
gases are expected to be largely dominated by the DDI. 

Deep ID optical lattices (formed by counter- 
propagating lasers) may slice a gas into non-overlapping 
samples. For non-dipolar (short-range interacting) par- 
ticles these non-overlapping gases may be considered as 
independent parallel experiments. The situation is com- 
pletely different in dipolar gases, since the DDI leads to 
inter-site interactions. For weak DDI (e.g. atomic mag- 
netic dipoles) these inter-site interactions lead to scatter- 
ing between BECs at different sites (lfj| : co llective exci- 



tations shared by non-overlapping sites II[ , or damping 
of Bloch oscillations || . For bosonic polar molecules the 
non-local dipolar effects may be much stronger, leading 
to fascinating effects as pair-superfluidity for ladder-like 
lattices [HI and filament BEC flj|. 

Filamentation is indeed an interesting possibility intro- 
duced by the DDI. This phenomenon, first suggested in 
the context of ferrofluids by de Gennes and Pincus [3], 
has attracted a considerable theoretical interest for the 
case of classical dipoles jlEj . Dipolar chains in classi- 
cal ferrofluids were recently observed in superparamag- 
netic iron colloids II6I1 and single-domain magnetite col- 
loids [13] ■ In Ref. [13] it was shown that a similar fila- 
mentation process may occur for bosonic polar molecules 



Figure 1: Polar fermionic molecules in a three- well poten- 
tial may remain unpaired, form fermionic trimers, or bosonic 
dimers between nearest neighbors or next-nearest neighbors. 



in deep lattices, which may organize into chains sustained 
by an attractive inter-site DDI, forming a dipolar chains 
liquid (DCL) which may Bose-condense 

In this Letter we consider DCLs of fermionic polar 
molecules. Far from being a trivial extension of the 
bosonic case, fermionic polar molecules lead to a very dif- 
ferent and rich physics. Whereas for bosonic molecules 
the chains are obviously bosons, for fermionic molecules 
the bosonic or fermionic character of the filaments de- 
pends on the number of molecules in the chain. This 
has particularly relevant consequences when the num- 
ber of available lattice sites is odd. Here we focus on 
the simplest non-trivial case, namely a three-well poten- 
tial (Fig. [I}. For simplicity we restrict our discussion 
to the ideal gas regime, where inter-filament interactions 
are neglected. Although this approximation is of lim- 
ited quantitative validity (and would demand mesoscopic 
samples in specific ID arrangements as discussed below), 
it contains already many of the qualitatively new features 
which may be expected for more general scenarios of po- 
lar Fermi molecules in deep ID and 2D optical lattices. In 
particular, the competition between trimer/dimer bind- 
ing and trimer Fermi energy results in a non-trivial de- 
pendence of the character of the DCL as a function of the 
number of molecules per site N. If N is smaller than a 
critical N c the DCL is a Fermi-degenerate gas of trimers. 
However for N > N c the trimers coexists with a Bose 
mixture formed by pseudo-spin- 1/2 dimers and spinless 
dimers, leading to a peculiar Bose-Fermi mixture. We 
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show that these Bose-Fermi DCLs may be probed by 
monitoring the spatial distribution of the molecules. 

We consider fermionic polar molecules with mass m 
and electric dipole d in a deep three-well potential along 
the z-direction, with inter-site spacing A. This arrange- 
ment may be created by e.g. selectively emptying all sites 
of a strong ID optical lattice except three neighboring 
ones. The potential barriers are large-enough to prevent 
any inter-site hopping. Due to reasons discussed below, 
the analysis of the problem simplifies notably if the gas is 
considered as strongly confined along the y-direction (e.g. 
in a single node of a lattice as that in the z-direction). 
Along the remaining x-direction we consider a shallow 
harmonic confinement with frequency lu. The molecules 
interact via the DDI Vd(r) = d 2 (l — 3 cos 2 0)/r 3 , where 
9 is the angle formed by r with the dipole orientation. 
We assume that the dipoles form an angle a with the 
z direction, such that sin 2 a = 1/3. Although this par- 
ticular orientation and the ID character of the sites are 
not needed for the formation of the DCL gas, which may 
occur also in stacks of 2D sites 13[, this particular sce- 
nario allows both for a strong attraction between dipoles 
placed on top of each other and for a vanishing DDI be- 
tween molecules at the same site. This largely reduces 
the inter-filament interaction, allowing for a simplified 
ideal gas scenario, as discussed below. 

The attraction between polar molecules placed on 
top of each other may be strong-enough to bind two 
or more polar molecules into self-assembled chains 
(Fig. Q]). Whereas for bosonic molecules these chains 



are in any case bosons 13(, for fermionic molecules the 
fermionic/bosonic character of the filaments depends on 
the odd/even number of molecules in a given chain. 
In particular, the three-well configuration allows for 
fermionic trimers (and of course monomers), and two 
different kinds of bosonic dimers, namely those between 
two molecules at nearest neighbors (type I dimers), and 
those between two molecules at the uppest and lowest 
site (type II dimers) (Fig. []}. Note that dimers I are 
actually pseudo-spin- 1/2 bosons, since dimers in sites 1 
and 2 are not equivalent to dimers in sites 2 and 3. 

The ground-state of a single filament of M molecules 
is calculated as for bosonic molecules [ljj, and hence we 
just sketch for completeness the basic ideas. Let r, (pj) 
be the position (momentum) operator of a molecule at 

E^r,/M, 



site j. Introducing P = X)-j=iPj> R 



Qj = Pj — P) s j = v j R = { x j, Vji z j}> the Hamiltonian 
splits into H = Hcm + Hreh where Hcm = P 2 /2Mm + 
Mmuj 2 R 2 /2 describes the filament center-of-mass and 
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approximated by a strong isotropic harmonic oscilla- 
tor of frequency lo±. The wavefunction of the j-th 



molecule is chosen as ipj(xj — Xjo)ipj(yj)(pj(z 



z jo), 



where fjirj) 
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exp ^— t)"/2Rq) j \J -\ZtvRq , with lj_ = y/h/mui_i and Rq 
is the variational width of the x wavepackets. For deep 
lattices one may approximate l± — > (energy corrections 
are < 1% for depths > 14 recoil energies h 2 ir 2 /2mA 2 ). 

Minimizing the energy of straight filaments (xjo = 0) 
with respect to Rq we obtain the filament binding energy 
[l8j |. We denote as — Et, —Er,,i and —Ez>,n the binding 
energies for, respectively, trimers, dimers I, and dimers 
II. These energies grow with the dipole strength Uq — 
md 2 /h 2 A. There exists a critical Uq such that for Uq < 
Uq the composites unbind (Rq > I ho = y/h/mui). Note 
that £/ *(T) < U$(D,I) < llf(D,II) (Fig. [2]) due to the 
different strength of the DDI in each composite. In the 
following we consider the regime Uo > Uq(D,II), where 
Rq <C Iho for all of the possible composites of Fig. [TJ 

Transverse filament excitations contribute to the gas 
entropy, being relevant at finite temperature T. In addi- 
tion, and contrary to the case of bosonic molecules [3] , 
transverse modes are important for fermionic molecules 
also at very low T since they may significantly reduce 
the trimer Fermi energy. For a chain of M molecules, we 
obtain the low-lying modes ^=1 m after expanding the 
chain energy E around its minimum, and diagonalizing 
d 2 E/dxjdxji , where j, j' = 1, . . . , M. 

In the following we consider the filament statistics, 
assuming an ideal filament gas. This largely simplifies 
the analysis of the problem, while allowing for the dis- 
cussion of key qualitative features of these systems, in 
particular the competition between different Bose and 
Fermi composites. This approximation is just quantita- 
tively valid for mesoscopic samples in the arrangement 
discussed above, as we discuss at the end of this Letter. 

The fermionic or bosonic character of the chains is re- 
flected by the average occupations for trimers, dimers I, 
dimers II and monomers: 
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where Ns,j denotes the average occupation of individual 
molecules at site j, £ 



VT;D,I;D,II 



the relative motion. The on-site yz confinement is 



the transverse filaments 
modes of the different composites, e„ = fku(n+ 1/2) the 
harmonic oscillator levels and (3 = 1/ksT the inverse 
temperature. In the previous expressions we have as- 
sumed symmetric configurations such that the number of 
dimers I in sites 1-2 is the same as the number of dimers 
I in sites 2-3, and equal to ND,i(n,fD,i)- Note that 
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Figure 2: Binding energy (in units of Eq = ft 2 /mA 2 ) of the 
different composites of Fig.[T]as a function of Uo = md 2 /ft 2 A. 

Ml = M3 is the chemical potential for molecules at the 
uppest and lowest sites, whereas (12 denotes the chemical 
potential for molecules in the middle site. These different 
chemical potentials are necessary to fulfill the normaliza- 
tion conditions, in which we assume N molecules per 
lattice site. Imposing symmetry between the uppest and 
the lowest sites, these conditions acquire the form: 

N = N T + N D ,i + N D>II + N s ,i, (6) 
N = N T + 2N DJ + N S .2, (7) 

where Nt, Nn,i, Nn,n, N$,i and Ns,2 denote respec- 
tively the total number of trios, dimers I in sites 1-2 (or 
2-3), dimers II, monomers in site 1 (or 3) and monomers 
in site 2. From © and ([7]) we obtain fii(N,T) and 
H2{N,T), and from (|2H5J) the occupation numbers. 

Due to the attractive DDI between molecules in the 
filament, the most bound chain is the trimer. The differ- 
ence in binding between dimers and trimers induces that 
for sufficiently small N and at low-enough T the DCL be- 
comes a degenerate Fermi gas of trimers. The trimers fill 
up oscillator levels (and also transverse trimer modes) up 
to the corresponding Fermi energy Ep(N), which equals 
Nhuj for rigid filaments but it is actually smaller due 
to the transverse trimer modes. However, if the num- 
ber of molecules per site is sufficiently large, the growth 
in Fermi energy overcomes the binding energy difference. 
This transition may be easily estimated by comparing the 
average energy per molecule for the case of two trimers 
and that for the case of 2 dimers I and one dimer II. This 
leads to a condition for the critical number of molecules 
per site N c (U ,u), E F {N C ) = 2E T - 3{E DJ + E DJI )/2 
(which we have confirmed numerically). Note that N c 
grows with growing Uo and decreasing uj. For N < N c 
the DCL is a degenerate trimer gas, whereas for N > N c 
the trimer gas coexists with a mixture of pseudo-spin- 1/2 
bosons (dimers I) and spin- less bosons (dimers II). 

The peculiar properties of the DCL translate into the 
spatial molecular distribution integrated over the three 
sites. For N < N c and N < £i t /gj, only trimers in their 
internal ground state are formed, and hence the gas be- 
haves as a spin-less Fermi gas of particles of mass 3m, pre- 
senting a Thomas- Fermi density profile (1 — (x/R) 2 ) 1 ^ 2 




Figure 3: Integrated density profiles of the molecules, for 
£i T /u < N < N c (top), and N > N c (bottom). We consider 
Uo = 2, lo/2-k = 1Hz, m = lOOamu, which lead to iV c = 230. 

with R/Iho = y/2N/3. For £_ 1t /uj < N < N c , the DCL 
is still a trimer gas, but transversal modes may be pop- 
ulated. In that case the density profile departs from the 
Thomas- Fermi profile (Fig. [3J top), due to the appear- 
ance of internally excited trimers in low harmonic oscil- 
lator levels. For N > N c the density profile changes dra- 
matically. Note that since we consider ID gases, dimer 
BEC is strictly speaking precluded. However, due to fi- 
nite size the dimers quasi- condense (at low-enough T) 
occupying the few lowest levels of the harmonic oscilla- 
tor. Hence when TV surpasses N c a Bose cloud nucleates 
at the trap center. As a result the distribution of the 
polar molecules shows a Gaussian-like peak at the trap 
center (Fig. 02 bottom). 

The DCL presents as well an intriguing finite temper- 
ature physics due to the role of filament modes and the 
different binding energy of dimers and trimers. This 
is particularly clear from a finite T analysis of a DCL 
with N < N c (Fig. |4|). Note that whereas at very low 
T the DCL is purely a trimer Fermi gas, at finite T it 
becomes more favorable to populate dimers than to pop- 
ulate higher excited trimer states. As a consequence the 
system presents a striking thermal enhancement of the 
bosonic modes. Interestingly, contrary to the standard 
situation, this leads to a maximal central peak density 
for a given finite T. For even larger T the central density 
decreases again due to the occupation of dimers at higher 
oscillator modes, and the breaking of the filaments into 
individual molecules. 

For N > N c and U > U£(D, II) the DCL is at low T 
a basically pure Bose gas of dimers I and II (except for a 
small trimer fraction). Since both dimers have mass 2m, 
the difference between them cannot be discerned from the 
analysis of the integrated density profile of the molecules. 
However the different binding energy and excited dimer 
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Figure 4: Temperature dependence of the fraction of 
molecules in trimers (solid), dimers I (dashed), dimers II (dot- 
ted) and monomers (dashed-dotted). We consider the param- 
eters of Fig. |3] with N < N c . 

modes for both types of dimers may be studied spectro- 
scopically to reveal the dual nature of the mixture. If 
N > N c but U£(D,I) < U < U£(D,II), dimers II 
are precluded, and hence the DCL will become at low 
T a Bose-Fermi mixture of dimer-I bosons and degener- 
ate monomers at sites 1 and 3 (which act as a pseudo- 
spin-1/2 fermions). Again, this exotic mixture could be 
revealed from the corresponding dual density profile. 

The discussed ideal gas analysis allows for a relatively 
easy understanding of key qualitative features charac- 
terizing fcrmionic polar molecules in deep optical lat- 
tices under more general conditions, as the competi- 
tion between filament-binding energy and Fermi energy 
of fermionic chains, the relevant role of the filament 
modes at zero and finite T, or the formation of pe- 
culiar mixtures of different types of composite bosons 
and fermions. However, the quantitative validity of the 
ideal gas approximation is rather limited (also for bosonic 
molecules [13|), even for the previously discussed ID ar- 
rangement with the particular choice for the dipole ori- 
entation. We may estimate the importance of the inter- 
filament interactions by comparing the inter-trimer in- 
teractions with the Fermi energy of rigid trimers (ep — 
Nfkv). For deeply bound chains (i?o ^ A/2) and at inter- 
filament distances x > A we may approximate the inter- 
action between molecules at different chains as that be- 
tween two point dipoles. Adding up these interactions we 
may estimate the mean inter-trimer DDI Vf / (x) , where x 
is the mean inter-trimer distance. For the case of d = 0.8 
Debye, m — 100 atomic mass units, A = 0.5/im, and 
lu/2tt = 1Hz, we obtain U ~ 2, and N c ~ 230. For this 
value x ~ 1.7A and V// ~ 0.33e_p. The ideal gas ap- 
proximation is hence quantitatively valid only for dilute 
mesoscopic samples (as those considered in our numerical 
calculations). Once the dimer Bose gas nucleates at the 
trap center the ideal gas condition is quickly violated, due 
to the larger bosonic densities, although the formation of 
the dual density profile (similar to that in Fig. [3|) still 
holds. For stacks of 2D sites the ideal gas approxima- 
tion fails even for extremely low 2D densities. However, 
the formation of dimer mixtures beyond a given critical 



density should also occur for 2D arrangements. These 
mixtures may be considered as weakly-interacting for 2D 
densities n such that nrl < 1 with r* = mcP/ti 2 . For 
the previous values this demands n < 1.1 X 10 8 cm~ 2 . 
For N > N c (and U > U^{D,II)) weakly-interacting 
dimers will form a BEC of three different bosons (dimers 
I in 1-2, 1 in 2-3, and II), whose properties will largely de- 
pend on the precise determination of the different inter- 
dimer interactions, which will be the subject of a future 
work. 

Summarizing, fermionic polar molecules in three-well 
potentials are expected to form a rather peculiar filament 
gas. Depending on the filling per site and the interac- 
tion strength we expect that the character of the chain 
gas ranges from a pure trimer gas at low fillings, to a 
bosonic mixture of pseudo-spin- 1/2 and spin-less dimers 
for large-enough fillings and dipole strengths. Note, fi- 
nally, that molecules in even larger number of sites may 
form a quantum gas mixture of increasing complexity. 
Dipolar chain liquids are hence an exciting perspective 
for on-going experiments with polar fermionic molecules. 
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